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ABSTRACT 

This effort extends into high frequency (>500 Hz), an 
earlier developed adaptive control algorithm for the 
suppression of thermo-acoustic instabilities in a liquid- 
fueled combustor. The earlier work covered the 
development of a controls algorithm for the suppression 
of a low frequency (-280 Hz) combustion instability 
based on simulations, with no hardware testing involved. 
The work described here includes changes to the 
simulation and controller design necessary to control the 
high frequency instability, augmentations to the control 
algorithm to improve its performance, and finally 
hardware testing and results with an experimental 
combustor rig developed for the high frequency case. 
The Adaptive Sliding Phasor Averaged Control 
(ASP AC) algorithm modulates the fuel flow in the 
combustor with a control phase that continuously slides 
back and forth within the phase region that reduces the 
amplitude of the instability. The results demonstrate the 
power of the method - that it can identify and suppress 
the instability even when the instability amplitude is 
buried in the noise of the combustor pressure. The 
successful testing of the ASPAC approach helped 
complete an important NASA milestone to demonstrate 
advanced technologies for low-emission combustors. 

INTRODUCTION 

Lean-burning, low emission combustors are being 
investigated for aircraft gas turbine engines. Lean 
combustion is shown to be advantageous for NOx 
emissions and efficiency, but suffers from an increased 
susceptibility to thermo-acoustic instabilities. These 
instabilities typically result from coupling of the 
fluctuating heat release 1 or coupling of the equivalence 
ratio modulation of the combustion process with the 
lightly damped acoustics of the combustion chamber. 


These instabilities cause pressure oscillations within the 
combustor that can reduce component life and potentially 
lead to premature mechanical failures' 

In recent years, suppression of thermo-acoustic 
instability has been attempted through active control. 2 
The goal of these active control efforts is to reduce the 
energy concentrated at the instability frequency and to 
reduce the overall amplitude of the combustor pressure 
oscillations. Some active control concepts involved 
speaker actuation or placing another heating 
source. 3,4,5,6,7,8,9 Speaker actuation simplifies the problem 
by eliminating the need for fuel modulation control with 
all its adverse process characteristics. However, acoustic 
actuation may not be directly applicable to aircraft-type 
engines due to the harsh environment and weight 
restrictions. Others have attempted to actuate the fuel 
modulation with some limited success. 10,11,12,13,14,15 These 
later control investigations have demonstrated some 
ability to reduce the frequency spectra of the instability 
but have been less successful at reducing the time 
domain pressure fluctuations. 

Many of these active combustion control (ACC) 
techniques rely on accurate knowledge of the plant 
dynamics. However, the combustion process is rather 
complicated and difficult to model. The combustion 
process characteristics include large dead-time phase 
shifts, large wideband noise compared to the amplitude 
of the instability, exponential growth of the instability, 
frequency and phase shift randomness, 16 and a system 
that transitions through inherently unstable operation. 
Because of these adverse process characteristics, the 
instability may grow very fast to unacceptable levels 
when even momentary loss of control tracking occurs. 
This momentary loss of control is inevitable when the 
instability has been suppressed down to the noise level. 
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This successful demonstration of ACC for the HF 
case was carried out with two distinct controllers 
developed by NASA. The second controller was 
observer based, called Multi-Scale Extended Kalman 
(MSEK) by Le 17 . Both controllers yielded good results. 
A summary of the whole NASA ACC effort and testing 
of the two controllers is reported by DeLaat. 18 

In the effort described in this paper, an adaptive phase 
shifting controller, 19 previously developed for 
suppression of a simulated Low-Frequency (LF) 275 Hz 
instability in a liquid-fueled single nozzle combustor, 
will be extended for applications in the High Frequency 
(HF) > 500 Hz combustion instability regime. The 
control algorithm senses pressure fluctuations due to the 
instability, filters, and phase shifts these pressure 
fluctuations to actuate the fuel modulation in order to 
generate combustion pressure or heat release oscillations 
that are opposing to the chamber acoustics. The control 
is based on adaptive phase shifting, where the phase of 
the control vector continuously bounces from one 
boundary of an effective stability region (i.e. a phase 
region that favors instability suppression) to the other, in 
a back and forth sliding motion. The Adaptive Sliding 
Phasor Averaged Control (ASP AC) algorithm requires 
very little detailed knowledge of the plant dynamics. 
Since this methodology is not model based, it can be 
expected that this controller can perform well at different 
operating conditions or at deteriorating engine 
conditions, where the structure of the engine dynamics 
can vary. 

The ASP AC controller is demonstrated on a liquid- 
fueled single-nozzle combustor rig. Previous work 
involving the development of the ASP AC approach will 
be extended here from the LF simulated configuration to 
the HF simulation and to control testing on the actual 
combustor rig. This paper documents changes made to 
the LF configuration in order to convert it to a HF 
simulation with corresponding control design changes, as 
well as augmentation of the controller to improve its 
performance and adaptability. The paper is organized as 
follows. A description of the combustor rig is given, 
followed by a description of the plant reduced-order 
model used for the HF control development. Changes or 
improvements to the control method required for control 
of the HF instability are then described. Some brief 
results of the simulation studies using the ASP AC 
method are presented. Combustor testing is discussed 
and results that demonstrate the effectiveness of the 
control method are presented. 

COMBUSTION INSTABILITY RIG 

In order to focus control development toward 
realistic combustion instabilities, a combustor rig which 



Figure 1 - Combustor Rig for Combustor 
Instability Control Research 

(Photo by United Technologies Research Center) 

replicates an engine combustion instability has been 
developed. 20,21 The magnitude of the pressure oscillations 
in the engine was sufficient to cause unacceptable 
vibratory stresses in the turbine. 

The rig successfully replicates the engine instability 
and operates at engine pressure and temperature 
conditions. The single-nozzle combustor rig has many of 
the complexities of a real engine combustor including an 
actual engine fuel nozzle and swirler, dilution cooling, 
and an effusion-cooled liner (Figure 1). For the 
conditions corresponding to the mid-power condition 
chosen for evaluation (T3=770 °F, P3=200psia, fuel-air 
ratio=0.03), test results established the existence of a 
combustion instability at approximately 566 Hz. 



Frequency (Hz) 

Figure 2 - Comparison of Engine and Baseline 
Combustor Rig Pressure Spectra 


NASA/TM— 2003-212535 


2 



A comparison between the pressure spectrum in the 
engine and in the single-nozzle combustor rig at 
comparable operating conditions is shown in Figure 2. 
As can be seen, the combustor rig approximates the 
engine instability frequency and amplitude. The 
combustor rig exhibits significantly higher wide band 
noise compared to the engine. As will be discussed later 
the noise level is the primary limiting factor of how 
much instability suppression can be achieved with a 
given control design. The level of noise in the rig makes 
it harder to extract the instability from among the noise 
and creates an additional challenge for the controller. 
The rig none-the-less provides a suitable, realistic test 
environment for combustion instability control research. 

The research combustor rig was developed in 
partnership with Pratt & Whitney and United 
Technologies Research Center (UTRC). Experimental 
testing with the combustor rig is taking place at UTRC. 

COMBUSTION INSTABILITY SIMULATION 


In a previous effort, a simulation for a LF (~ 275 Hz) 
instability, observed with a different rig configuration, 
was developed and simulations and control development 
were carried out. 19 Because the instability frequency of 
the HF rig configuration is more engine like, the decision 
was made to adjust the controller design in order to 
attempt control of the HF instability. The basic ASP AC 
control/simulation structure remains unchanged for 
control of the HF instability, (Figure 3). Shown is the 
self-sustained oscillations loop ( G , , G A , NL) and the 
controller with the fuel valve actuation. In this closed 
loop simulation the overall combustor pressure is 
expressed as the sum of the instability pressure generated 
by the self-excited flame and acoustics loop, and the 


pressure caused by the controller actuated fuel 
modulation. There are many variations in previous 
published works as to how the physics of this control 
problem are interpreted and simulated. However, this 
controlled process boils down to a vector summation, 
vector rotation problem. 19 

The gain and damping for the reduced-order model 
of the LF and the HF combustor simulations remain the 
same. The frequencies were changed to simulate the HF 
instability. Based on this change, the new flame 
dynamics and acoustics models are as follows: 


G f = 


4500 

s + 3770 


( 1 ) 


G 


A 


-3.26e6 

s 2 +510.84s + 1.305e7 


( 2 ) 


Based on these choices for G F and G A the instability is 
self-excited, that is, it requires no input other than non- 
zero initial conditions (via fuel injection modulation). 
The G a transfer function dynamics correspond to an 
instability frequency of 575 Hz with low damping 
(g = 0.0707 )• The instability loop response time is ~ 0.6 
sec based on instability response measurements. 2 The 
sampling time of 10 kHz was selected to reflect the 
sampling frequency of the test hardware, which defines 
the controlled phase resolution. Given the choices of 
damping and sampling time this simulation is 
numerically unstable (limited by the saturation non- 
linearity NL), which means the open-loop response also 
depends on the sampling time. Numerical instability is 



Filter 

Figure 3 - Combustion Instability Control Block Diagram 
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Figure 4 - Combustion Rig Filtered Instability 
Pressure (psi) versus Reference Sine Wave 


not an undesirable feature with this simulation, as long as 
the results are repeatable (which they are). The saturation 
non-linearity and the combined model (G v ) of the fuel 

valve dynamics, feed lines, injector, equivalence ratio 
modulation, and flame dynamics models are as follows: 


NL = tanh( p) 

(3) 

25.3e6 

(4) 

s 2 +6e3^ + 25.3e6 


where, p is the acoustic pressure. The purpose of the NL 
is to provide a soft limit for the instability. By soft limit 
the meaning here is that the amplitude of the instability 
will be limited when no forcing is applied, but also the 
amplitude can grow when forcing or fuel modulation is 
applied. Without NL, and due to the low damping and the 
numerical instability, the amplitude of the simulated 
combustor instability will grow unbounded. For the 
transfer function q no modeling information was 

available, therefore, the assumption was that it had a 
second order, somewhat damped response and a pass- 
band that covered the instability frequency of ~575 Hz. 
Based on that, G v has the same damping and gain as for 

the LF configuration, only its frequency was changed. 

The noise level was adjusted based on the averaged 
wideband noise of the rig in Figure 2, which produced 
averaged noise ~ 7 to 10 times higher than the instability. 
The open loop plant simulations produced similar results 
to those described in Reference 19 except that in this 
(HF) case, the time domain response is dominated by 
noise. This is due to the high noise-to-signal ratio of the 
HF rig configuration compared to the LF configuration. 

As before, two identical band pass filters are placed in 
series in the feedback path. These filters achieve two 
objectives: 1) filtering out the process and sensor noise; 


and 2) representing the dead-time phase shift in the 
combustion process. Normally, the time delay associated 
with the delivery, injection, atomization, vaporization, 
and burning of the fuel would be included in the dead- 
time phase shift of the plant. However, for convenience, 
this dead-time phase shift has been included in the 
feedback path. Note : From block diagram this is not a 
standard feedback control system. The filter structure 
and attenuation characteristics for the HF case are the 
same as in Reference 19, only the band pass frequency 
was changed. Based on that, the state space model of the 
HF band pass filter used in the simulation, with a pass 
band frequency of 510 to 640 Hz, is as follows: 

x = Ax + Bu 


With 



" -0.8 

0 

0 

3.58 

0 

0 


0.8 

-0.8 

-0.8 

0 

3.58 

0 


0 

0.8 

0 

0 

0 

3.58 

,4 = 1000 

-3.58 

0 

0 

0 

0 

0 


0 

-3.58 


0 

0 

0 


0 

0 

-3.58 

0 

0 

0 

B = [800 0 

0 0 

0 of 

c = [o 0 

1 0 

0 o] 


This filter produces the same phase shift at the instability 
frequency as before. However, phase shift associated 
with dead time delay is expected to double with the HF 
rig configuration (assuming line lengths remain the 
same), because the frequency is approximately doubled. 
Therefore, with the HF configuration the total phase shift 
could be in the order of 1400 degrees. 

Examining combustor rig experimental data shows 
the HF combustor instability exhibits net phase walks 
that can be as severe as 400 degrees per control cycle of 
40 Hz (control cycle will be explained in the next 
section). In addition, this instability frequency is random, 
varying within approximately +/-10 Hz, centered at 
approximately 545 Hz. Also, the instability continuously 
varies in amplitude. Figure 4 shows the rig instability 
pressure, filtered to remove the noise, against a reference 
sine wave of the same frequency and constant amplitude. 
This figure shows the instability undergoing amplitude 
modulations as well as random phase walk. With a more 
coherent instability, where the noise-to-signal ratio is 
low (like in the case of the LF rig configuration used 
previously or the actual engine — Figure 2), these adverse 
affects are less pronounced. Based on this, it is expected 
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that suppressing a coherent instability with low noise-to- 
signal ratio would be a lot easier to accomplish. Analysis 
and testing to confirm this are being carried out and 
expected to be published at a later date. 

The ASP AC was also evaluated briefly against a more 
physics-based model of the combustion rig. 22 The 
physics based model was programmed in FORTRAN 
while the controller was in C language. It is expected that 
more evaluations of this type will be conducted in the 
future. 


CONTROL ALGORITHM 


Examining the control structure of Figure 3 a little 
more carefully gives a better understanding of one of the 
major challenges facing this control design. As the 
instability pressure is being suppressed via the controlled 
fuel modulation, the amplitude of the required fuel 
modulation continues to decrease. At some point the 
modulation amplitude becomes comparable to the 
amplitude of the noise. Therefore, as the instability is 
being suppressed, noise increasingly drives the process 
which ultimately will cause loss of control phase 
tracking (phase adaptation) resulting in the process 
becoming inherently unstable. Inherently unstable in this 
sense means that at some point as the instability is being 
suppressed, there will be no controlled phase modulation 
that will cause additional suppression. Rather, any 
modulation at this point, no matter what its amplitude 
and phase, will in fact cause amplification of the 
instability. Therefore, the desired operating command in 
this feedback control system, which calls for zero steady- 
state instability amplitude, is an inherently unstable 
point. In that sense, noise with this process becomes the 
limiting factor of how far down the instability can be 
suppressed. Of course, other process characteristics such 
as large phase delays are important too, but none are as 
important as noise in this control design. 

Although the combustion process is rather complex, 
from a controls perspective, the desire is to simply find 
an opposing modulation phase that will suppress the 
instability. From this viewpoint, this behavior is 
equivalent to vector summation and vector rotation. The 
ASP AC algorithm senses the instability pressure, filters 
it, and averages the signal over a period of time (i.e. 
reduces the control bandwidth) in order to overcome the 
large process phase delay. It then phase shifts this signal 
to control the fuel valve in order to suppress the 
instability. The ASP AC controller consists of the 
following parts: 

1. Continuous phase shift (10 kHz sample time) 

2. Adaptive control phase shift (40 Hz) 



Figure 5 - Bode Plot of Band-Pass Filter 

3. Discontinuous exponential gain modulation control 
(-500 Hz) 

4. Adaptation of control parameters ( 4 Hz) 

Continuous Phase Shift (part 1) 

In part 1 , the instability is sampled at a rate of 1 0kHz 
and filtered to remove the noise. Then a phase shift angle 
(as computed in part 2) is applied, and this signal is 
output as a command to the fuel actuator, also at 10kHz, 
to control the fuel modulation. When the controller is 
used to control the plant instead of the simulation, the 
design of the band pass filter is changed to reduce the 
phase shift. This is because the filter in the simulation (as 
discussed before), was also used to simulate the plant 
phase shift, which is no longer needed for the actual rig. 
The state space model of the band pass filter design for 
the actual HF plant controller, with a pass band of 450 to 
700 Hz is: 


x- Ax + Bu 

y = Cx 


With 


.4 = 1000 


-2.2627 

1.6 

-3.51 

0 


- 1.6 

0 

0 

-3.51 


3.51 

0 

0 

0 


0 

3.51 

0 

0 


( 6 ) 


B = [l600 0 0 Of C = [0 1 0 0] 


The Bode plot of the band pass filter is shown in Figure 
5. When the controller was used against the simulation, 
the controller performed better as the width of the pass 
band was increased. However, as will be discussed in the 
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controller testing section, the filter bandwidth had to be 
adjusted for the rig. 

Adaptive Control Phase Shift (part 2) 

This part of the control calculates the phase angle of 
the filtered instability signal (used in part 1), necessary to 
suppress the instability. Therefore, in part 2 the same 
filtered instability signal sensed at a 10 kHz rate is 
averaged over a 40 Hz cycle in order to determine 
whether the amplitude of the instability signal is 
decreasing or increasing. If the instability amplitude is 
decreasing, the controller steps and advances the phase in 
the same direction as before. If the signal amplitude is 
increasing, the controller reverses the direction it is 
advancing the phase angle. 

This part the algorithm, as shown in Figure 6, also 
computes a restricted control region. The control activity 
is restricted to operate within this region. The motivation 
for computing this restricted control region is to know in 
advance the phase region that promotes instability 
suppression. The reason for this is because 
intermediately in this control process, control phase 
tracking will be lost as suppression approaches the wide- 
band noise level. A new restricted control region is 
computed whenever control suppression has been lost for 
a certain successive number of control cycles, and then 
suppression subsequently is recovered and continues for 
a predetermined number of control cycles. 

Within the restricted control region the controller 
operates, sliding back and forth within a so-called 
effective stability region. The boundaries of this stability 
region are defined as the last phase angle in either 
direction that causes instability suppression. Once the 
controller phase stepping exceeds the phase angle at the 
boundary of this stability region, causing instability 
amplification, the controller reverses direction. This 


stable region continues to shrink as the instability is 
being suppressed further, until there exists no phase 
angle that will cause further suppression. This becomes 
the case because at this point the instability feedback 
signal is small enough that the noise now becomes the 
dominate driver in the process. When this inherently 
unstable point is reached the instability begins to 
amplify, and eventually the controller again succeeds in 
reestablishing control, thereby repeating the process over 
and over. 

Discontinuous Exponential Gain Modulation Control 
(part 3) 

This feature is active only when a restricted control 
region has been established (part 2). In this portion of the 
algorithm the controller proportional gain, by which the 
sensed instability is multiplied in part 1 , is modulated on 
and off. Even though, the on/off frequency is not fixed, 
the smallest duration of time for the gain “on” condition 
is equal to the pass-band filter low-pass corner 
frequency. The condition that triggers the gain “on” 
condition is decreasing instability amplitude over this 
time duration. The gain “off’ condition is triggered by an 
increase in the instability signal amplitude. Instead, 
triggering the gain “on” when the instability amplitude is 
increasing and “off’ otherwise, also worked well. This 
tends to indicate that generally discontinuous control 
modulation may be better than continuous. When the 
decision is made to turn the gain “on”, a small time delay 
is applied first, and then the gain decays exponentially 
starting from an initial value. A simulation of the 
controller proportional gain is shown in Figure 7. This 
initial value of gain is also adapted in part 4 of the 
algorithm. 

What motivates this Discontinuous Exponential Gain 
Modulation Control (DEGMC), is the very nature of a 
resonant mode coupled with noise — that is, when the 


Pressure from 
Fuel modulation 


Boundary of effective 
Stability region 



Overall combustor pressure 


Pressure from 
instability 


Boundary of restricted 
control region 


Figure 6 - ASP AC Phasor Diagram in a Stationary Frame of Reference 
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CONTROLLER TESTING 


1 .3 1 .32 


Figure 7 - Controller Proportional Gain Showing 
Discontinuous Exponential Gain Modulation 


instability is sufficiently suppressed to the point that the 
process noise, Figure 3, provides a significant part of the 
excitation, additional suppression may become futile 
(reaching the point where the system becomes inherently 
unstable). The unintended result is amplification of the 
instability. Suppressing the instability only while its 
amplitude is decreasing seems feasible considering that 
resonant modes undergo amplitude modulation and any 
amplification will not be monotonically increasing on a 
cycle-by-cycle basis. 

The motivation for the exponential gain decay is to 
overcome in a sense the large phase delay and its adverse 
effects that cause amplification of the effective control 
gain. The phase delay, which can be on the order of 1000 
degrees for the HF configuration, can cause the effective 
control gain to vary by nearly an order of magnitude 
from its lowest to its highest value. 

Adaptation of Control Parameters (part 4) 

In part 4, the control parameters are adapted 
sequentially and continuously at a rate of 4 Hz or less. 
The original controller design for the LF configuration 
discussed in Reference 19 did not include this control 
parameter adaptation. The control parameters which are 
adapted are the controller proportional gain, the 
exponential gain and the exponential gain delay 
discussed in part 3, and the frequency of the adaptive 
phase shift discussed in part 2. In the case of the 
proportional control gain, this part of the control 
algorithm adapts the initial value of the gain used for 
exponential decay, (Figure 7). The adaptation 
methodology used in this case is called Adaptive Binary 
Search. 23 This adaptation methodology essentially steps 
the control parameter adjustment in the direction that 
reduces the averaged instability signal amplitude, and 
reverses direction, halving the step size in the process 
when the amplitude of the signal increases. 


The objective of the combustor rig testing was to 
assess the overall performance of the controller, in its 
ability to suppress the HF engine-like instability. Two 
separate tests were carried out using the ASP AC to 
control the HF configuration of the combustor rig 
described in section 1 . The plan was to test each of the 
controller parts discussed in the previous section, and to 
tune the various controller parameters. However, due to 
facility problems during both tests, the test time was 
significantly limited. Also, some of the parameters were 
not successfully recorded, thus the controller did not get 
the chance to be exercised thoroughly in a methodical 
way. Despite these limitations, the data obtained 
successfully validated the control concept and 
demonstrated the ability of the controller to suppress the 
HF instability. The controller worked well even though 
the rig instability was not very coherent and there was 
considerable combustion noise. 

The order of the controller testing was as follows: 
test the ability of the controller to suppress the instability 
and find the band pass filter that produces best 
performance, exercise the DEGMC portion of the 
algorithm, and finally exercise the parameter adaptation 
part. 

Initial Instability Suppression and Filter Testing 

At the beginning of rig testing a baseline case was 
run without feedback controls to establish the 
uncontrolled instability level. This baseline instability 
case was repeated after each series of controller testing. 
For the initial closed-loop control testing, only part 1 and 
part 2 of the ASP AC were run with band-pass filters of 
differ bandwidth. The filters tested in this portion of the 



Figure 8 - Filtered Combustor Rig Pressure versus 
Filtered and Phase Shifted Fuel Valve Control 
Command 
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Frequency (Hz) 

Figure 9 - Instability Suppression with and 
without Discontinuous Exponential Gain 
Modulation Control 

controller evaluation had the following pass bands; no 
filter, 350 to 800, 400 to 750, 450 to 700, 510 to 630 and 
540 to 604 Hz. The pass band of the last filter tested was 
slightly outside the instability center frequency (the 
instability frequency was not precisely known before the 
test). Through this part of testing, it was established that 
the band pass filter discussed in the previous section of 
450 to 700 Hz pass band, achieved the best performance 
in suppressing the instability. The assumption is that with 
a relatively narrower band pass filter (narrower than 450 
to 700 Hz), the noise has a more pronounced effect since 
it’s been concentrated near the instability frequency. On 
the other hand, too wide a pass band will not perform as 
well (as can be seen in Figure 2, the rig has high noise 
power at lower frequencies). There is also an added 
benefit with the fact that the wider pass band of 450 to 
700 Hz is better for robustness: That is, this wider band 
is expected to cover the variability of the instability 
frequency for all operating conditions of a given engine 
configuration. Therefore, with the use of such a filter, 
there is no need to extract the instability frequency in 
real time. As seen in Figure 2, the shape of the actual 
engine spectral density differs from that of the rig, which 
means that the best filter design for the actual engine 
could be somewhat different from the simulation or the 
rig. 

Figure 8 shows the combustor rig pressure 
and the controlled fuel modulation command. Both 
signals are filtered to show only the instability and the 
respective control command. Because the process has 
large phase delay, the control command has also been 
processed with a fixed phase shift. If the control signal 

is represented as u c - A sin(a>t + 6 ) , then the processed 
control signal is u c p roc - A sin {cot + 6 + </>) , where 


<p = Rem {</> p dela y / 2n) and <t> p4day is the plant 

phase delay in radians. Rem stands for remainder. With 
this preprocessing, the combustor pressure and the 
control command shown in Figure 8 are displayed as if 
the plant has no phase transport delay, and the opposing 
phase of the controlled modulation to the combustor 
pressure is evident. This confirms that there is actual 
control of the instability taking place, as opposed to the 
controller causing some sort of excitation of the 
combustor pressure with the indirect affect of a net 
reduction in the instability pressure. This figure also 
shows that when the instability pressure is suppressed to 
nearly zero, the phase of the combustor pressure changes 
drastically, driven by the noise which dominates at this 
point. Even though the controlled signal phase starts 
opposing combustion pressure soon after the combustor 
pressure begins to grow again, the pressure still 
continues to grow for while. This sort of behavior 
supports earlier discussions that the plant becomes 
inherently unstable at suppression levels near the noise 
level. Second, the instability growth persists even though 
the controller is continuously opposing this action. This 
behavior seems to go beyond the effects of the large 
dead-time phase delay of the plant, and rather it points 
out that in addition, there exists instability inertia or 
some self-reinforcing mechanism. This later result is 
preliminary and will be further investigated in the future. 

Discontinuous Exponential Gain Modulation Control 
Testing 

DEGMC was exercised during control testing on the 
HF rig. Figure 9 shows instability suppression on the rig 
with and without DEGMC. The additional suppression 
due to DEGMC is >10%. However, when put into 
perspective of the noise effects, as will be discussed later 
in this section, the significance of this additional 
suppression will become clear. Since this new control 
design methodology (DEGMC) has not been fully 
evaluated, more performance assessments will need to be 
done. 

After filter and DEGMC tuning, the amplitude 
spectra of the uncontrolled combustion instability against 
the controlled instability from the two separate tests are 
shown in Figure 10 and Figure 11. The amplitude 
suppression in the instability for test 1 and test 2 is 
approximately 30% and 40%, respectively. The low 
frequency oscillation seen in the first test, Figure 10, is 
due to a separate valve mean flow controller. This was 
fixed in the second test. 

From Figure 12, it can be seen that the combustion 
pressure is dominated by random noise. Because of the 
large amount of noise, the controller could only reduce 
the time domain RMS pressure by a few percent. But, 
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Figure 10 - Test 1, Amplitude Spectral Density of 
Uncontrolled versus controlled instability 



Figure 11 - Test 2, Amplitude Spectral Density of 
Uncontrolled versus controlled instability 


even if complete elimination of the instability were 
possible, it would make only a few percent reduction in 
the amplitude of the time domain signal. Since the noise 
floor of the actual engine combustion pressure is 
significantly less than that exhibited by the rig, Figure 2, 
it is expected that in the case of the actual engine the 
controller will be able to more significantly reduce the 
peak time domain amplitudes and RMS value of the 
combustion pressure. 

Since noise is the primary factor limiting the amount 
of instability suppression that is achievable, a new 
measure is proposed here for comparing the effectiveness 
of different controller designs in suppressing combustion 
instabilities. Instead of comparisons involving 
percentage or dB reduction compared to the uncontrolled 
instability amplitude, this preferred measure eliminates 
comparisons involving the uncontrolled instability. 
Rather, in terms of amplitude spectral density, the 
measure of success in suppressing combustion instability 
is based on how close to the noise floor the controller can 
suppress the instability. This new measure for comparing 
the effectiveness of various controller designs, defined 
below, is the Instability Suppression Ratio (ISR): 


S, 

Where S N is the average amplitude of the wide band 
noise spectra in the region of the instability, and S/ is the 
peak amplitude spectra of the suppressed instability 
(provided that the controlled instability is discernibly less 
than the uncontrolled case). The interpretation of the 
average noise amplitude S N , in cases where the wide 
band noise is not flat, would be the noise value at which 
the instability would be indistinguishable from the noise 
floor. Based on this definition, an ISR 1 of 1 would be the 
maximum possible suppression achievable, in which case 
the instability will be indistinguishable from the noise on 


a spectral density plot. Based on this measure, the values 
of ISR 1 achieved during the first and second test are 
approximately 0.70 and 0.60, respectively. The fuel 
valve that was found to be partially clogged during the 
second test plus the limited test time could be reasons for 
the reduced performance of the controller during the 
second test. As discussed before and under this new 
measure of success, even 10% improvement in the ISR 1 
is significant, especially when the ISR 1 approaches 1, 
making it increasingly more difficult to achieve 
additional suppression. 

Adaptation of Control Parameters Testing 

Control parameter adaptation was exercised during 
the second test. However, because of the limited time 
available and because there were problems with the data 
recording, an evaluation of this part of the control 
algorithm could not be conducted. Assuming that the 
instability signal is averaged for sufficiently long time, 
and that the adaptation steps are sufficiently large to 
avoid possible drifts of the optimized parameters, 
parameter adaptation is expected to work. 



Figure 12 - Combustion Pressure of Controlled 
Instability 
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FUTURE PLANS 

Future plans include more testing with the HF rig 
configuration for a more thorough testing of the control 
algorithm. This high frequency testing will hopefully 
allow a more methodical tuning of the control parameters 
for better performance. Also, this testing will allow a 
more thorough evaluation of the parameter adaptation 
portion of the algorithm. Furthermore, frequency sweep 
tests will provide data to help derive a more accurate 
process model for future controls development. 

Future plans also include testing the ASP AC 
algorithm against the extended, LF rig configuration. The 
LF configuration uses the same combustor with added 
length and volume in the diffuser section. The LF 
configuration exhibits a strong coherent instability at 
approximately 280 Hz. This test will help evaluate the 
ability of the control algorithm to adapt to changing 
conditions, preferably without the need for major code 
modifications. This result has already been shown, at 
least partially, when the controller, developed against a 
LF simulation, was used without any essential code 
modification to control the HF combustor rig 
configuration. Also, since prior control developments 
have been done against the LF configuration, the LF test 
will allow a direct comparison with earlier control 
methods work. 24 

It is expected that future experimental work will be 
preceded by testing the ASPAC against the physics 
based model. 22 The objective is to use the physics based 
model for validation purposes before hardware testing. 

Finally combustion instability test data for the LF 
and HF rig configurations as well as actual engine data 
are being analyzed to assess the characteristics of these 
different combustor configurations to gain a better 
understanding of instability behavior. 

CONCLUSION 

A combustion instability suppression control 
algorithm has been developed and successfully 
demonstrated on a combustor rig that incorporates many 
of the complexities of a real engine. The Adaptive 
Sliding Phasor Averaged Control (ASPAC) algorithm 
performs well with relatively high combustion noise, 
high dead-time phase shift, net random phase walk, some 
frequency randomness and amplitude modulation of the 
instability waveform. In addition, the algorithm responds 
relatively fast to suppress the instability and also to 
recover from loss of control tracking. This loss of 
tracking is unavoidable when the instability is 


sufficiently suppressed during control action. The control 
algorithm combines tracking control of the relative phase 
between instability pressure and fuel modulation with 
phase adaptation and control parameter adaptation. This 
combination is designed to allow good control of the 
instability even in cases when there are changes in the 
phase required for instability suppression such as those 
due to net random phase walks. Also, since this 
controller is not model based, its applicability is expected 
to be wide in terms of operating conditions and engine 
degradation, and even directly applicable to engines with 
different or variable geometries. Future work includes 
demonstration of this control approach on a different 
combustor rig configuration, integration of the control 
algorithm with a physics based model, and development 
of emission reduction controls combined with 

combustion instability controls. 
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